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0.  Introduction 

In  recent  years,  there  has  been  considerable  interest  in 
determining  the  manner  in  which  complicated  oscillatory  phenomena 
can  occur  in  dynamical  systems  through  bifurcations  (see,  for 
example,  [5] , [6] ,  [9] , [10] , [11] , [12] , [18] ) .  In  this  paper,  we 
consider  the  same  type  of  problems  for  low  dimensional  systems. 

More  specifically,  in  Section  2,  we  suppose  f  is  a  vector 
field  on  1R"  which  corresponds  to  a  bifurcation  point  of  degree 
one;  that  is,  the  vector  field  possesses  either  a  saddle-node  or 
a  degenerate  focus  or  a  degenerate  periodic  orbit  or  a  homoclinic 
orbit.  The  first  problem  that  is  considered  is  the  characterization 
of  the  flow  for  a  perturbation  of  this  vector  field  when  the  per¬ 
turbation  is  allowed  to  be  periodic  in  the  independent  variable. 

By  using  the  c.lassical  theory  of  integral  manifolds  appropriately, 
we  are  able  to  give  a  complete  solution  for  the  saddle-node  or 
degenerate  focus.  The  theory  also  applies  to  the  degenerate  focus 
even  when  the  perturbations  are  almost  periodic.  For  the  cases 
corresponding  to  a  degenerate  periodic. orbit  or  a  homoclinic  orbit, 
only  partial  results  are  obtained. 

In  Section  3,  we  consider  the  behavior  of  the  solutions  for 
an  autonomous  equation  in  dimension  >  2  near  an  equilibrium 
point  when  the  linear  variational  equation  has  more  than  two  eigen¬ 
values  on  the  imaginary  axis.  More  specifically,  we  consider  an 
equation  in  with  two  purely  imaginary  and  one  zero  characteris¬ 

tic  root.  Using  the  theory  of  normal  forms  and  polar  coordinates. 
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one  can  consider  this  as  the  periodic  perturbation  of  an  autonomous 
equation  in  IR  .  We  first  give  a  complete  bifurcation  diagram  for 
the  autonomous  equation  under  generic  hypotheses  on  the  quadratic 
and  cubic  terms.  Each  point  on  a  bifurcation  curve  is  a  bifurca¬ 
tion  point  of  degree  one.  The  theory  of  Section  2  is  then  applied 
to  obtain  a  complete  description  of  the  problem  in  1R^.  This 
generalizes  results  of  [5], [9], 

We  give  explicit  approximations  for  the  bifurcation  curve 
and  show  the  specific  structure  of  the  bifurcation  near  the  homo¬ 
clinic  orbit.  Near  the  Hopf  bifurcation,  our  results  could  also 
be  obtained  from  known  results  on  behavior  of  diffeomorphisms 
(see,  for  example,  [15],  [19]).  However,  our  methods  based  on  the 
more  classical  theory  of  integral  manifolds  permits  the  considera¬ 
tion  of  perturbations  which  are  even  almost  periodic. 

•  In  Section  3,  we  also  discuss  an  equation  in  1R*  with  two 

pairs  of  purely  imaginary  roots.  Using  the  theory  of  normal  forms 

and  polar  coordinates,  the  problem  is  reduced  to  the  discussion 

2 

of  an  autonomous  equation  in  1R  perturbed  by  higher  order  terms 
involving  two  angle  variables.  This  corresponds  to  a  perturbation 
of  the  equation  in  IR^  into  ]R  .  The  autonomous  equation  in  1R 
involves  no  quadratic  terms.  Under  some  generic  hypotheses  on  the 
terms  up  through  fifth  degree,  we  give  a  bifurcation  diagram  for 
the  autonomous  problem. which  has  the  property  that  each  bifurcation 
is  of  degree  1. 

In  particular,  we  prove  the  conjecture  in  [9]  that  there  is 
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a  first  integral  of  the  approximate  system  at  the  critical  values 
of  the  parameters  where  a  Hopf  bifurcation  is  conceivable.  This 
has  also  been  recently  obtained  independently  in  [6],  From  the 
consideration  of  higher  order  terms,  we  resolve  the  nature  of  the 
bifurcation  at  this  point.  Again,  using  the  analysis  of  Section 
2,  we  mention  partial  results  for  the  bifurcations  in  1R  . 

1 .  Summary  of  known  results 

If  ft  is  an  open  set  in  IR^,  r  *  3ft,  and  *  ft  U  r, 
let  <%2  be  the  set  of  Cr-vector  fields  from  to  1R^ 

which  are  transverse  to  T.  Two  vector  fields  f,g  are  equivalent 
f  ^  g,  if  there  is  a  homemorphism  on  ST  which  maps  orbits  of 
f  onto  orbits  of  g  preserving  the  sense  of  time.  An  f 
is  structurally  stable  if  there  is  a  neighborhood  U  of  f  such 
that  g^f  if  g  €  U.  An  f  is  a  bifurcation  point  of 

degree  one  if  f  is  not  structurally  stable  and  there  is  a  neigh¬ 
borhood  U  of  f  such  that  g  6  U  implies  g  is  either  struc¬ 
turally  stable  or  g  f. 

Theorem  1.1  (see  [2], [17]). 

An  f  €#2  is  structurally  stable  if  and  only  if  every 
equilibrium  point  and  every  periodic  orbit  is  hyperbolic  and  there 
are  no  connections  between  saddle  points. 


The  set  of  bifurcation  points  of  degree  one  also  is  complete¬ 
ly  classified.  The  basic  result  is  the  following  (see  [1] , [20] , [21]) . 
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Theorem  1.2.  A  vector  field  f  £•£?!;,  r  —  is  a  bifurcation 
point  of  degree  1  if  and  only  if  there  is  a  neighborhood  W  of 
f  and  a  submanifold  r  of  codimension  one  in  W  such  that 
W  ^  T  «  Uj  U  U2  where  each  g  €  Uj  is  structurally  stable  but 

g^h  if  g  £  Uj,  h  €  For  g  £  r,  only  one  of  the  following 

situations  prevails: 

(i)  g  €  r  has  an  elementary  saddle-node  at  Xp,  there  are 
no  equilibrium  points  of  g  near  Xp  if  g  £  and  a  saddle 
and  node  near  xQ  if  g  €  U2« 

(ii)  g  6  T  has  an  elementary  focus  at  xQ,  there  is  no 

periodic  orbit  of  g  near  xQ  iff  g  €  U,  and  a  periodic  orbit 

near  xQ  if  g  £  U2  -  the  generic  Hopf  bifurcation. 

(iii)  g  £  T  has  a  periodic  orbit  y  which  is  stable  from 
one  side,  unstable  from  the  other,  g  €  has  no  periodic  orbit 
neat'  y  and  g  €  U2  has  two  hyperbolic  periodic  orbits  near  y. 

(iv)  o0  «  tr  3f(0)/3x  f  0,  g  £  T  has  a  homoclinic  orbit 
containing  a  saddle  point  Xp,  g  £  has  a  saddle  near  xQ  and 
no  periodic  orbit  near  y,  g  £  has  a  saddle  point  and  a  unique 
hyperbolic  periodic  orbit  near  y  which  coalesce  as  g  r. 

Each  of  the  cases  (i)-(iv)  are  shown  in  Figure  1. 
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2.  Periodic  perturbations. 

For  vector  fields  of  dimension  >  2,  the  literature  on 
structural  stability  and  bifurcation  is  extensive,  but  the  results 
are  less  complete.  We  illustrate  some  of  the  difficulties  by 
considering  perturbation  of  a  planar  autonomous  system  by  periodic 
terms. 

Suppose  X , y  are  real  parameters  and  consider  the  system 

(2.1)x>u  x  *  f (x,X)  +  g(t,x,vi) 

where  g(t,x,y)  =  g(t+l,x,y) ,g(t,x,0)  =  0. 

We  assume  that  f(*,0)  is  a  bifurcation  point  of  degree 

one.  More  precisely,  suppose  U  =  Uj  ij  r  'j  U2  as  in  Theorem  1.2, 

f(-,0)  €  T,  fC-,X)  €  Uj  for  X  <  0,  f(*,X)  €  U2  for  X  >  0  and 

the  curve  (f(-,X),X  e  1R)  is  transverse  to  r  at  X  =  0. 

The  vector  field  f(*,X)  is  structurally  stable  for  X  f  0 

and  the  manner  in  which  the  change  is  made  from  one  type  of 

structurally  stable  system  to  another  at  X  *  0  is  one  of  the 

situations  shown  in  Fig.  1  and  described  in  Theorem  1.2. 

By  an  application  of  the  Implicit  Function  Theorem,  for  any 

Xq  €  1R,  there  is  an  n(Xg)  >  0,  v(Xq)  >  0  such  that  each  hyperbolic 

equilibrium  point  of  the  autonomous  equation  (2.1).  n  is  perturbed 

A0’ 

to  a  1-periodic  solution  of  (2 . 1 ) x  y  for  |X-XQ|  <  nCwg)* 

|y|  <  v(Xq).  The  classical  theory  of  integral  manifolds  (see,  for 
example,  [7])  implies  any  hyperbolic  periodic  orbit  is  perturbed  to 
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a  smooth  hyperbolic  invariant  manifold  in  1R  *  1R  which  is  like 
a  cylinder  whose  cross  section  is  one  periodic  for  |X-XQ|  <  n(XQ), 
|y|  <  v(Xq).  Identifying  the  cross  section  at  t  =  0  with  the 
one  at  t  =  1,  we  obtain  a  torus.  Thus,  if  XQ  f  0;  that  is 
f(*,X0)  is  structurally  stable,  then  the  flow  for  the  perturbed 
equation  is  understood.  Of  course,  the  size  of  the  perturbation, 

| X - X Q |  <  n(XQ),  |u|  <  v(Xq),  for  which  this  is  true  depends  upon 
Xq .  The  difficulty  arises,  therefore,  in  a  neighborhood  of 
Xq  =  0.  Let  us  discuss  this  situation  in  more  detail. 

Suppose  case  (i)  occurs  for  X  =  0,  y  =  0.  Let  p  be  the 
saddle  node  in  Fig.  1  and  let  v  be  a  unit  vector  in  the  null 
space  of  3f(p,0)/3x.  Using  the  classical  method  of  Liapunov- 
Schmidt,  one  obtains  a  bifurcation  function  G(a,X,y)  for 
|a|,|X|,|y|  sufficiently  small  so  that  (2.1),  has  a  1-periodic 

a  j  y 

solution  x(t)  in  a  sufficiently  small  neighborhood  of  p  if  and 

only  if  there  is  an  (a,X,y)  such  that  G(a,X,y)  =  0  and  the 

constant  a  is  such  that  x(t)dt  =  p  +  av.  Furthermore,  G(a,0,0) 
22  J° 

aa*  +  0(|aj  )  as  |a|  0,  where  a  ■  is  a  nonzero  constant.  Using 

this  fact,  it  is  not  too  difficult  to  show  that  a  neighborhood  U 
of  (X,y)  =  (0,0)  can  be  decomposed  as  U  =  U^  'J  T  U  Uj  where  T 
is  a  smooth  curve  containing  (0,0)  such  that  the  period  one  map 
of  the  Eq.  (2.1)..  in  a  neighborhood  of  p  is  the  same  as  the 

A  >  V 

one  in  Fig.  1.1.  This  gives  a  complete  solution  of  the  problem 
in  case  (i) . 

Suppose  f(*,0)  satisfies  (ii)  of  Theorem  1.2;  that  is,  the 
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where  G(t,x,X,y)  is  1-periodic  in  t. 


G(t,x,X,y)  =  F(<KX,v)  (t)+x,X)  +  g  Ct ,  CX  »p)  (t)+x,y) 

(2.6) 


-  F(<j>(X,y)  (t)  ,X)  -  g(t,<J>(X,y)  (t)  ,y)  • 


Thus,  G(t,0,X»y)  =  0  for  all  t  €  1R,  (X,y)  €  V,  G(t,x,X,0)  = 
F(x,X). 

One  now  can  introduce  polar  coordinates  x  ^  (p,0), 
x^  *  p  cos  0,  x2  =  -P  sin  0,  to  obtain  the  equations 


(2.7) 


X,V 


0  =  ojx  +  0(t,0,p,X,y) 
p  =  axp  +  R(t , 0 , p , X ,y)  . 


Eq.  (2.7).  must  be  analyzed  for  (p,X,y)  in  a  neighborhood 
of  zero.  To  use  the  theory  of  normal  forms  or  averaging  effectively, 
we  make  the  further  nonresonance  hypothesis 


(2.8)  mw0  +  n/2ir  f  0,  m,n  integers,  |m|  +  |n|  £  4. 

Hypothesis  (2.8)  implies  there  is  a  change  of  variables 

0  9  +  u(t,9,p,X,y)  ,  p  h-  p  +  v(t,0,p,X,y) 


with  u,v  periodic  in  t  and  6  respectively  of  period  1  and 
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2ir/u)Q,  u  being  0(|p|),  v  being  o(|p|),  for  X,y  =  0, 
such  that  Eq.  (2.7),  becomes 

A  ,  y 

e  =  0JX  +  8(X,y)  +  0(t,8,P,X,y) 

(2-9h,M 

p  =  (ax+Y(X,y))p  +  d(X,y)p3  +  R(t,0,P,X,y) 

where  0  vanishes  for  p  =  0,  R  =  o(|p|3)  as  P  -*  0.  Further¬ 
more,  the  properties  of  G  in  (2.6)  and  the  fact  that  f ( * , 0) 
corresponds  to  a  generic  Hopf  bifurcation  implies 

(2.10)  6(X,0)  =  0,  Y(X , 0)  =  0,  d(0,0)  t  0,  aQd(0,0)  <  0, 

the  latter  inequality  being  the  case  if  the  picture  is  the  one 
shown  in  Fig.  1  for  case  (ii). 

Let  C  be  the  curve  defined  in  a  neighborhood  of  (X,y)  =  0 
by 

(2.11)  C  =  { (X,y) :  ax  +  y(X jy)  =  0}  . 

This  curve  is  well  defined  since  f  0  and  Y(X,0)  =  0. 
Introduce  new  small  parameters  (e,6),  e  =  ax  +  y(X,y),  6  =  y  so 
that  C  corresponds  to  e  »  0. 

To  be  specific  in  what  follows,  suppose  that  d(0,0)  <  0. 

In  this  case,  for  (X,y)  e  C,  that  is,  t  =  0,  the  integral 

manifold  p  ■  0  of  (2.9),  is  uniformly  asymptotically  stable 

A ,  y 


11 


for  5  sufficiently  small.  For  e  <  0,  the  same  is  true  and 
one  can  actually  find  an  eg  >  0  and  a  6g  such  that  this  is 
true  for  -eQ  <  c  <  0,  |6|  <_  5Q.  For  each  e  >  0,  the  theory 
of  integral  manifolds  shows  there  is  an  invariant  torus  which 
is  uniformly  asymptotically  stable  for  |5|  <  6g(e),  and  attracts 
everything  in  a  neighborhood  of  p  =  0  except  the  line  p  =  0. 

This  fact,  together  with  the  stability  at  e  *  0  implies  we  can 
find  positive  Cg,6g  such  that  all  of  the  above  assertions  are 
valid  for  |e|  <  Cg,  J  5  J  <  6Q.  Details  of  the  proof  are  not  given 
here.  We  note  that  by  introducing  new  variables  p  =  6p  and 
e  *  6c,  we  may  apply  the  integral  manifold  theorem  in  [7] 
directly  to  obtain  the  desired  tori.  We  may  also  prove  the 
uniqueness  of  such  tori  by  using  the  uniqueness  method  in  [4]. 

These  results  are  summarized  in  the  following  theorem. 

Theorem  2.1.  Suppose  f(*,0)  satisfies  (ii)  of  Theorem  1.2  and 
the  equilibrium  point  of  f(*,X)  is  stable  for  A  <  0  and 
conditions  (2.4)  and  (2.8)  are  satisfied.  Then  there  is  a  neigh¬ 
borhood  W  o_f  x  =  0,  a  neighborhood  V  of  (A,y)  *  0  and  a 
smooth  curve  C  c  V  containing  (0,0)  such  that  V  *  u  C  U  V2 

with  Vj ,V 2  open  connected  sets  such  that  the  following  conclusions 
hold  for  all  solutions  x(t)  of  (2.1).*  ..  with  initial  values  in 

- -  -  ,  A,y  - 

W: 

(i)  For  (A,p)  6  Vj,  there  is  a  1-periodic  solution  <|>(A,y) 

of  (2.1).,  ..  with  both  characteristic  multiples  less  than  one  and 
—  a  ,  y  ' 

x(t)  -  4>(A,y)(t)  -*-0  as  t  +  «. 


(ii)  For  ( X ,  y)  €  C,  there  is  a  1-periodic  solution 
4>(X,y)  such  that  both  characteristic  multipliers  are  one  and 
the  same  conclusion  as  in  (i)  holds . 

(iii)  For  (X,y)  £  there  is  a  1-periodic  solution 

<J> (X , p)  of  (2.1).,  ..  with  both  characteristic  multipliers  greater 

a  f  y 

than  one  and  an  invariant  cylinder  in  1R^  with  1-periodic  cross 
section  which^uniformly  asymptotically  stable  and  x(t) ,  x(0)  f 
<J>(X,y)(0),  approaches  this  cylinder  as  t  -*■  ®. 

If  the  equilibrium  point  of  f(*,0)  is  unstable  and  f(*,0) 
satisfies  (ii)  of  Theorem  1.2,  then  a  similar  result  is  obtained 
with  obvious  changes  in  the  stability  statements. 

Remark.  The  above  theorem  is  also  true  for  quasi-periodic  g 
provided  the  appropriate  conditions  on  the  frequencies  in  terms 
of  roots  of  unity  are  assumed.  Similar  remarks  apply  when  0  is 
a  vector  and  all  functions  are  periodic  in  each  component  of  e. 


Now,  suppose  that  f(’,0)  satisfies  (iii)  of  Theorem  1.2. 
This  case  is  much  more  complicated  than  the  previous  one  because 
there  is  no  way  to  obtain  any  solution  whose  qualitative  properties 
are  known  for  all  X,y  (in  the  previous  case,  we  used  strongly 
the  fact  that  a  1-periodic  solution  was  known  to  exist  for  all 
X,y).  On  the  other  hand,  some  results  can  be  obtained,  but  they 
are  not  as  complete  as  in  case  (ii). 


For  X  =  y 


has  a  unique  periodic  orbit 


“  0,  Eq.  (2 . 1) Q  ^  q 
Yq  which  is  asymptotically  stable  from  one  side  and  unstable 
from  the  other.  Suppose  Yq  ■  {x0(9),  9  €  1R)  where  Xg(t)  is 
a  periodic  solution  of  (2.1)0  Q  of  least  period  2ir/a>0,  u>0  >  0. 

If  p  represents  distance  from  Yq  along  a  normal,  then  one 
can  introduce  a  local  coordinate  system  around  Yg  which  takes 
x  h-  (e,p)  to  obtain  the  equations 


1  +  0(t,0,p,X,y) 


(2.2) 


X,y 


p  *  R(t,0,p,X,y) 


where  0,R  vanish  for  (p,X,y)  *  0,0,0),  are  1-periodic  in  t 
and  '  2ir/oi0  periodic  in  8. 

Let  us  suppose  that  we  can  average  to  obtain  an  equivalent 
set  of  equations 

0  *  1  +  aX  +  by  +  0(t,6,p,X,y) 
p  =»  CX  +  dy  +  ep  +  fp^  +  R(t,9,p,X,y) 


C2-13)x,y 


where  a,b,c,d,e,f  are  constants,  0,R  for  p  =  0  are 
0((|x|  +  |y  |)2)  as  X  ,y  -*>  0  and  R  for  X  *  y  -  0  is  0(p3). 

The  hypothesis  that  f(*,0)  satisfies  (iii)  of  Theorem  1.2  implies 
cf  <  0,  e  ■  0.  Since  c  |>  0,  let  us  introduce  new  parameters  (6,e) 
as  6  •  y,  e  *  cX  +  dy  to  obtain  the  equations 
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0  •  1  ♦  a'6  ♦  b'e  ♦  0’ (t,e,p,e,6) 

(214>e,{ 

p  *  t  ♦  fp2  ♦  R' (t,0,p,e,6) 

with  0',R'  having  the  same  type  of  order  relations  as  before. 

For  any  e  small  enough  such  that  cf  >  0,  it  is  shown 

there  is  a  6g(e)  >  0  such  that  no  solution  of  (2.14)g  remains 

in  a  neighborhood  of  p  *  0  for  all  t  if  0  <  6  <  6Q(e).  This 

is  the  analogue  of  X  <  0  and  the  picture  in  Fig.  1  for  the  period 

one  map.  If  ef  <  0,  then  the  theory  of  integral  manifolds  shows 

there  is  a  6Q(e)  >  0  such  that  there  are  two  invariant  cylinders 

1  2 

T£  Te  of  (2.14)£  fi,  0  <  6  <  50(e),  with  periodic  cross 

9  9  12 
section;  that  is,  two  invariant  tori  T  . ,  T  x.  One  of  these 

c ,  o  e  ,  o 

tori  is  completely  unstable  and  the  other  is  uniformly  asymptoti¬ 
cally  stable.  Furthermore,  any  solution  of  (2.14)  which  remains 

in  a  sufficiently  small  neighborhood  of  p  *  0  for  t  >  0  (t  <  0) 

1  2 

must  approach  either  T  .  or  T‘  .  .  as  t  -*•  »(t  -<*>).  Thus,  we 

G  f  0  G  y  0  < 

have  complete  knowledge  of  the  solutions  in  a  subset  of  a  neighbor¬ 
hood  of  (e,6)  ■  (0,0).  The  subset  is  depicted  in  Fig.  2. 


6 

I 


I 

I 

Figure  2. 


As  e  0,  the  number  6g(e)  may  a*so  approach  zero. 

1  2 

The  reason  for  this  fact  is  that  the  invariant  tori  T  .,  T  t 

C  |  O  E  ,  o 

for  a  fixed  e  may  fail  to  exist  as  smooth  invariant  sets  for 

some  small  j.  This  is  a  general  fact  in  differential  equations. 

Given  a  smooth  hyperbolic  invariant  torus  for  a  differential 

equation,  the  existence  of  a  smooth  invariant  torus  for  a  small 

perturbation  of  this  equation  generally  requires  that  the  strength 

of  the  hyperbolicity  is  greater  than  the  rate  of  attraction  of 

orbits  on  the  torus.  If  this  is  not  the  case,  the  invariant  sets 

develop  cusps  at  rational  rotation  numbers.  For  the  sets  T* 

T2  .  above,  the  strength  of  the  hyperbolicity  approaches  zero  as 
C  |  o 

c  0  for  any  fixed  6.  Thus,  it  is  conceivable  that  the  sets 

touch  for  some  e  without  coalescing  in  a  uniform  manner.  This 

is  WJby  we  cannot  complete  the  picture  in  Fig.  2.  The  behavior  in 

the  unshaded  region  will  be  extremely  complicated  and  no  one  has 

made  any  progress  toward  an  explanation.  The  classical  theory  of 

integral  manifolds  gives  no  information. 

Suppose  now  that  f(*,0)  satisfies  condition  (iv)  of 

Theorem  1.2.  This  situation  has  been  discussed  much  more  in  the 

literature  because  so  many  new  phenomena  occur.  If  p  is  the 

saddle  point  of  the  homoclinic  orbit  y  of  (2.1)q  q,  then  the 

Implicit  Function  Theorem  gives  a  1-periodic  solution  $(X,y)  of 

def 

(2.1)>u  for  X,u  small,  *(0,0)  -  p,  and  PXy  »  ♦(X,y)(0) 

will  be  a  saddle  point  for  the  period  one  map  it  •  irXy.  Let 
W*  ,  ,  wV  ..  be  the  stable  and  unstable  manifolds  of  P,  as  a 

*  f M  A , V 

fixed  point  of  *. 
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A  point  q  is  homoclinic  to  P.  if  q  E  W?  n  wV  . 

A  ,  H  *(|j  A,|J 

It  is  transverse  homoclinic  to  P.  if  it  is  homoclinic  to 

a  » y 

P.  and,  in  addition,  W*  ,  wV  are  transversal  at  q.  For 
Afy  a f y  a  y  y 

X  *  p  ■  0,  any  point  q  €  y  is  homoclinic  to  pQ  ■  xQ,  but  it 
is  not  transverse  homoclinic.  Under  a  small  perturbation  in  an 
appropriate  direction,  one  would  expect  that  W®  and  w!f 
intersect  transversally  at  some  point  q  close  to  y.  The 
behavior  of  the  solutions  of  (2.1).  is  then  very  complicated 
near  y.  In  fact,  one  can  show  there  are  infinitely  many  dis¬ 
tinct  periodic  solutions.  Also,  there  is  an  invariant  set  of 
v 

n  for  some  integer  k  for  which  the  flow  is  equivalent  to  the 
shift  automorphism  on  two  symbols  (see  114]). 

It  has  also  been  shown  that,  in  a  neighborhood  of  the 
bifurcation  point,  there  are  generically  in  g  infinitely  many 
hyperbolic  sources  (or  sinks)  (see  (16]).  This  suggests  that 
the  bifurcation  to  transverse  homoclinic  points  might  be  related 
to  some  other  types  of  successive  bifurcations  of  periodic  orbits 
of  longer  and  longer  period.  The  existence  of  these  successive 
bifurcations  have  not  been  studied  wheri  f(>,0)  satisfies  (iv)  in 
Theorem  1.2. 

For  the  case  where  f(’,0)  is  Hamiltonian,  successive  bi¬ 
furcations  appear  to  be  the  general  rule  as  indicated  by  the 
results  in  [3].  Let  us  summarize  these  results  for  a  special  case. 
Consider  the  equation 
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(2*15)a,u 


x  ■  y, 

y  -  x  -  x2  -  \y  *  yf(t) 


where  X,y  are  small  parameters,  f  is  continuous  and  f(t*l)  « 
f(t).  Eq.  (2.15)q  q  has  a  homoclinic  orbit  y*  but  it  does  not 
satisfy  (iv)  of  Theorem  1.2  since  the  divergence  of  the  vector 
field  at  the  saddle  point  (0,0)  is  zero.  Let 
Y  «  ((0,0)}  u  ( (p(t) ,p(t)) ,  —  <  t  <  -}  where  (p(t),p(t))  is 
a  solution  of  (2.15)0  Q  which  approaches  zero  as  t  -►  ±<*>.  We  may 
normalize  p  so  that  p(0)  *  0. 

Let 


(2.16)  h($)  -  j  p(t)f(t-*)dt/J  p2(t)dt 

-OO  -00 

and  suppose  that  every  extreme  point  of  h  is  a  strict  local 
maximimum  or  local  minimum.  Let  £  (0,1)  correspond  res¬ 

pectively  to  the  absolute  maximum  and  minimum  of  h.  We  now  state 
the  following  result  from  (3).  Parts  (i),  (ii)  were  also  essen¬ 
tially  given  in  [13]. 

Theorem  2.2.  Under  the  above  hypotheses,  there  is  a  neighborhood 

U  of  y  in  1R2,  a  neighborhood  V  of.  (X,p)  -  (0,0)  in  1R2, 

k  k 

an  integer  Y.q  >  0  and  curves  C^,  Cm  ill  V,  k  ■  kQ,kg+l , . . .  ,• 
(including  »)  which  have  parametric  representations  with  twice 


!*• 
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k  k 

continuously  differentiable  functions,  each  C^,  Cffl  contains 

(0,0),  cjj  -►  cJJ,  C*  -  C*  as  k  ♦  -  uniformly,  CjJ[C"]  is 

tangent  to  the  curve  defined  by  X  *  h($M)y  (X  *  h($m)y]  at 

(X,y)  ■  (0,0)  and  for  each  integer  k  €  [kg,*],  the  set 

V  V  (cjjj  U  C*)  »  Sj  U  Sj,  Sj,S2  open,  disjoint,  and  the  following 

conclusions  hold  for  Eq.  (2.15), 

_______________________  a  ,  y 

(i)  X  €  S“  implies  no  homoclinic  point  exists  in  U 

OQ 

(ii)  X  €  S2  implies  a  transverse  homoclinic  point  exists 

•  -OO 

—  s2 

(iii)  I_f  k  <  ®,  X  €  S  implies  no  subharmonic  solution  in  U 

(iv)  I_f  k  <  «,  X  €  implies  there  are  at  least  2k 

hyperbolic  subharmonics  in  U. 

k 

(v)  I_f  X  €  f\>|{  S2,  then  there  are  infinitely  many 

hyperbolic  nodes  or  foci  in  U  corresponding  to  fixed  points  of 
iterates  of  the  period  one  map. 

The  above  theorem  shows,  in  particular,  the  following.  If 
(X,y)  is  a  point  of  bifurcation  to  transverse  homoclinic  orbits; 
that  is,  (X,y)  €  or  C^,  and  W  is  any  neighborhood  of  this 
point,  then  there  are  infinitely  many  subharmonic  bifurcations  in 

w. 


Let  us  outline  the  proof.  In  (2.15).  ,  if  x(t)  ■  p(x+$)  ♦ 

*  t  M 

z  (x+$) ,  t  ♦  $  *  t,  then 


(2.17) 


z  ♦  a(t)z  -  F(t,z  ,z  ,X  ,y  ,♦) 
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where  a(t)  *  -1  +  2p(t),  the  derivative  of  -x  ♦  x2  evaluated 
at  p(t).  The  function  F  is  given  by 

F  =  -Az  -  Ap  +  uf (t - d>)  -  z2. 

The  first  step  of  the  proof  is  to  obtain  conditions  for 
the  existence  of  a  homoclinic  orbit.  This  is  equivalent  to  find¬ 
ing  a  solution  z  of  (2.17)  which  approaches  zero  as  t  -*•  ±». 

We,  therefore,  need  an  analogue  of  the  Fredholm  alternative  for 
bounded  solutions  on  1R  of  the  nonhomogeneous  equation 

z  +  a(t)z  =  F(t) 

where  F(t)  is  bounded  and  continuous  on  JR.  A  careful  analysis 
of  the  variation  of  constants  formula  shows  that  such  a  solution 
exists  if  and  only  if 


p(t)F(t)dt  »  0> 

•  oo 


If  this  condition  is  satisfied,  then  the  solution  is  unique. 

After  having  this  Fredholm  alternative,  one  can  adapt  the 
reduction  principle  of  Liapunov-Schmidt  or  the  method  of  alternative 
problems  to  this  situation  to  obtain  a  bifurcation  function 
G(<fr,A,y)  defined  for  $  €  1R,  (A,y)  in  a  neighborhood  of  (0,0), 
which  has  the  property  that  there  is  a  homoclinic  orbit  in  a  small 
neighborhood  of  y,  if  and  only  if  G($,A,y)  *  0.  Furthermore, 
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it  is  transverse  homoclinic  if  and  only  if  3G(<j>, X, y)/ 3<j>  t  0. 

To  solve  this  equation,  one  observes  that 

G(<M,y)  =  -X  4  yh(<fr)  +  0((|X|4|y|)2) 

as  X,y  -*•  o  where  h ( <t>)  is  defined  in  (2.16).  The  proof  of 
parts  (i)  and  (ii)  in  the  theorem  are  now  supplied  without  much 
difficulty. 

To  prove  the  other  parts,  one  uses  the  results  in  [8]  for 
obtaining  the  curves  of  bifurcation  to  subharmonics  of  order  k 
coming  from  the  orbit  of  period  k  of  (2.15)0  Q.  One  then  shows 
by  a  nontrivial  argument  that  the  conditions  on  h(a)  imply  the 
ones  in  [8]  uniformly  in  k.  This  will  complete  the  proof  (see 
[3]  ‘for  details) . 


3.  Two  purely  imaginary  and  one  zero  root. 

In  this  section,  we  consider  a  differential  equation  in 
1R^  for  which  the  linear  variational  equation  near  the  equilibrium 
point  zero  has  two  purely  imaginary  eigenvalues  and  one  zero  eigen¬ 
value.  To  simplify  the  situation,  it  will  be  assumed  also  that  a 
certain  type  of  symmetry  prevails.  More  specifically,  consider 
the  equation 


x  «  A(X)x  4  f (x,y) 


(3.1) 


y  -  By  ♦  g(x,y) 


2 


where  A, 8  are  small  real  parameters, 
4 

C  -functions, 


x  6  1R2 ,  y  €  1R,  f  ,g 


are 


A(A1  - 


(3.2) 


f(x,y)  »  0(|x|(|x|+|y|)) 


g(x,y)  =  0((|x|+|y|)2) 


as  I x  | ,  |y|  -*•  0.  The  hypothesis  on  f  implies  that  the  symmetry 
condition 


(3.3)  f (0,y)  «  0 

is  satisfied. 

Since  f(0,y)  =0,  it  is  legitimate  to  introduce  polar 
coordinates  for  x  *  (Xj,x2)  as  x^  =  r  cos  8,  x2  =  -r  sin  0. 

If  this  is  done  and  t  is  replaced  by.,  8,  one  obtains  the  equations 


(3.4) 


p  3  Ap  +  R(0,p,y,A,8) 

y  *  By  +  Y(e,p>y,A,6) 


where  R  *  0(| p| (| p| +|y| )) ,  Y  -  0( ( | p | + | y | ) 2)  as  p,y  ♦  0. 

The  problem  is  to  determine  the  behavior  of  the  solutions 
of  (3.4)  in  a  neighborhood  of  (p,y)  •  (0,0)  for  (A, 8)  in  a 


JL 
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neighborhood  of  (0,0).  To  discuss  (3.4),  one  can  use  the  theory 
of  normal  forms  (or  averaging)  to  transform  (3.4)  to  a  system 


(3.5) 


P  =  p(X+ay)  ♦  0 ( | p | (p2+y2) ) 
y  =  By  *  by2  +  cp2  +  0((|p|+|y|)3) 


as  p,y  -*•  0.  The  generic  situation  is  to  have  a,b,c  nonzero 
constants . 

To  simplify  the  computations,  and  also  to  consider  the 
most  interesting  case,  we  suppose  a=2,b=c=-l  so  that 
the  equations  are 


(3.6) 


p  =  p (A+2y)  +  0( ( p [ (p2+y2) ) 
y  -  By  -  y2  -  p2  +  0(( | p | + | y | ) 3) 


as  p,y  0. 

It  is  convenient  to  introduce  the  change  of  variables 
y  *-*■  8/2  +  y,  A  +  8  *-*■  a,  to  obtain  the  more  symmetric  form 


(3.7) 


p  =  P (a+2y)  +  0((p| (p2+(B/2+y)2) 

2 

Y  *  V  -  y2  -  P2  +  0((|p|+B/2+y)3) 


as  p,8,y  0. 
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If  we  perform  the  scalings 

(3.8)  p  -*•  ep,  y  ey,  8  ■+•  e,  a  -*■  ea,  t  -*■  e_1t 


the  new  equations  become 


(3.9) 


P  =  2py  +  ap  +  0  C | e | ) 
y  =  \  -  y2  -  P2  +  Odd J 


as  e  -*■  0. 

Equations  (3.9)  must  be  discussed  for  all  p  0 ,  y  e  1R, 
a  €  1R  and  e  in  a  neighborhood  of  e  =  0. 

For  a  =  0,  e  *  0,  the  function 


(3.10) 


V(p,y)  =  J-p 


is  a  first  integral  of  (3.9).  Since  the  Jacobian  of  the  vector 
field  in  (3.9)  has  trace  a,  it  follows  that  there  can  be  no 
periodic  orbit  of  (3.9)  unless  a  =  0.  This  remark  makes  much 
of  the  discussion  of  (3.9)  very  simple  for  a  i  0,  e  =  0.  In 
fact,  the  topological  structure  of  the  flow  is  determined  by  the 
qualitative  properties  of  the  equilibrium  points. 

The  equilibrium  points  of  (3.9)  for  e  =  0  are  p  *  0, 

1  2  2 
y  *  ±  j  for  all  a  and  the  point  y  *  *a/2,  p  *  (1-a  )/4  for 
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ct2  <  1.  For  a2  t  1,  the  points  p  *  0,  y  *  ±  1/2  are  hyper¬ 
bolic  saddles  or  nodes.  In  fact,  for  a  >  -1  (a  <  -1),  the  point 
(1/2,0)  is  a  saddle  (stable  node).  For  a  >  1(<1),  the  point 
(-  1/2,0)  is  a  saddle  (unstable  node).  For  a  >  0  (<0) ,  the 
..■point  (-a/2,  (l-a2)^^2/2)  is  a  stable  (unstable)  focus .  For  a  =  1, 
the  point  (-  1/2,0)  is  a  saddle-node.  For  a  =  -1,  the  point 
(1/2,0)  is  a  saddle-node.  These  facts  imply  that  the  phase 
portrait  for  (3.9)  for  a  /  0,  e  =  0  are  the  ones  shown  in  Fig.  3. 
The  only  bifurcation  points  are  when  a  =  ±1  and  these  are  of 
saddle -node  type.  For  a  =  0,  e  =  0,  the  function  V  in  (3.10) 
is  a  first  integral  and  the  phase  portrait  is  determined  from  the 
level  curves  of  this  function  shown  in  Fig.  3. 

As  remarked  in  the  previous  section,  for  each  a  /  0,  there 
is  an  e0(a)  -such  that  the  same  portraits  are  valid  for  the  period 
map  (remember  the  terms  in  e  are  periodic  in  0)  of  the  complete 
equations  (3.9)  for  0  1  e  <  eQ(a).  This  gives  3  uniform  estimate 
in  e  as  long  as  a  remains  in  a  compact  set.  To  obtain  a 
uniform  estimate  for  a  large,  we  return  to  the  original  equation 
(3.7)  and  introduce  the  scaling 

pn-ep,  yn-ey,  a^-e,  0  h-  e$ ,  t  *-►  e  _1t 

to  obtain  the  equations 


P  =  2py  +  P  +  0(|  e  | ) 
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In  this  equation,  we  need  to  discuss  only  the  case  where  e  is 
near  zero  and  6  varies  in  a  neighborhood  of  zero  because  the 
case  for  large  6  is  the  same  as  the  previous  one  with  the 
scaling  (3.8).  For  e  ■  0,  0  *  0,  the  equations  have  only  one 
equilibrium  point  p  =  0,  y  *  0  and  this  is  a  saddle -node.  The 
analysis  in  the  previous  section  applies  to  this  case. 

This  shows  that  the  only  case  remaining  to  discuss  in 
Eq.  (3.9)  is  a  =  0,  e  *  0.  The  manner  in  which  the  phase 
portrait  in  Fig.  3  changes  as  we  vary  (a,e)  in  a  neighborhood 
of  (0,0)  cannot  be  determined  with  specific  knowledge  of  the 
terms  in  order  e.  Thus  suppose  system  (3.6)  is  written  as 

p  *  p(A+2y+dp2+ey2)  ♦  0( ( | p | + |y | )4) 

(3.11)  - 

;•  y  =  By  -  y2  -  p2  +  fp2y  +  gy3  +  0((|p|+|y|)4) 


as  p,y  0. 

The  final  result  depends  on  all’.of  the  constants  d,e,f,g. 
To  simplify  the  computations,  we  suppose  d=l,  e=f*g*0. 

The  general  case  is  treated  in  the  same  way  although  more  computa¬ 
tions  are  involved.  Thus,  our  system  is 

P  “  p (A+2y+p2)  +  0((|p|  +  |y|)4) 

y  -  By  *  y2  -  p2  +  0((|p|+|y|)4). 


(3.12) 
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Making  the  same  transformations  as  before  and  using  the  scaling 
(3.8),  one  obtains  the  Equation  (3.9)  in  the  more  explicit  form 


(3.13) 


p  *  2py  +  ap  +  ep3  +  0(|e|2) 

y  -  \  -  y2  -  p2  +  0(|c|2) 


as  e  -*•  0. 

2 

Our  first  objective  is  to  analyze  Eq.  (3.13)  when  the  e 
terms  are  neglected;  that  is,  the  equation 


(3.14) 

»■- 


2py  +  ap  +  ep4 


For  a  =  e  =  0,  the  phase  portrait  is  given  in  Fig.  3.  Let 
(p0  (t)  »y0  (t) ) ,-  PQ  (t)  >  0,  be  the  solution  of  (3.14)  for 
a  *  e  *  0  describing  the  heteroclinic-orbit ;  that  is,  pg(t)  +■  0 
as  t  -*•  ±°o,  yQ(t)  -  1/2  as  t  -*•  ®,  y0(t)  1/2  as  t  -*•-«.  We 

may  assume  yQ (0 )  =0.  In  a  neighborhood  of  a  ■  e  *  0,  there 
will  be  a  curve  with  the  property  that  Eq.  (3.14)  has  a  hetero¬ 
clinic  orbit  joining  (-  1/2,0),  (1/2,0)  for  each  (a,e)  belong¬ 
ing  to  this  curve.  The  determination  of  this  curve  follows  along 
the  same  ideas  as  in  the  case  (iv)  of  the  previous  section  as  an 
application  of  the  method  of  Liapunov-Schmidt .  One  constructs  a 
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scalar  function  G(a, c) , . G(0 ,0)  ■  0,  such  that  there  is  a 
heteroclinic  orbit  if  and  only  if  G(a,e)  »  0.  Furthermore, 

G(a,e)  -  a  |  pg(t)dt  +  e  j  pj(t)dt  +  0(  (  ja|  +  |  e  | ) 2) 

as  a,e  -*•  0.  The  Implicit  Function  Theorem  implies  there  is  a 
unique  solution  a  -  a*(e)  of  G(a,e)  *  0  for  e  sufficiently 
small,  a*(0)  *  0  and 


6 


dgf  da* (0)_ 
3e 


r  pjctjdt 

J  .00 

f  P«(t)dt 
J  -  00 


1 

T  • 


The  Specific  value  of  6  as  -  ^  is  obtained  in  the  following 
way.  Since 


p0dy0  *  ^0  *  p0y0  p0)dt 

■L.c»dt  ■  iL°°dt 


0  *  PoV(pO»y(P  ■  *  p0y0  '  3^0 


we  have 
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from  dp0  =  2p0y0dt  and  the  symmetry  in  the  equation.  Since 
V(P0(t),y0(t)}  *  0,  we  have  P2(0)  «  3/4,  PZQ  *  3(l/4-y^), 
yQ  *  (l/4-pg/3)1^2.  Using  these  relations  we  have,  for  e  *  Pq 


Thus,  6  *  -1/2, 

The  function  a  *  a*(e)  satisfying  G(ct*(e),e)  *  0  for 
|ej  <  Cq  has  the  property  that  there  is  a  heteroclinic  orbit  of 
Eq.  (2.14)  for  (a,e)  *  (a*(e),e). 

Near  the  value  a  ■  0,  e  -  0,  '.there  is  also  a  curve  where 
a  Hopf  bifurcation  occurs.  For  a  ■  0,  e  *  0,  Eq.  (3.14)  has 
an  equilibrium  point  y  ■  0,  p  -  1/2.  The  Implicit  Function 
Theorem  implies  there  is  an  equilibrium  close  to  this  one  for 
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(a,e)  small  and  it  is  given  approximately  by  yQ  ■  -(a+e/4)/2, 
pQ  *  1/2  +  yQ.  Analyzing  the  stability  properties  of  this 
solution,  one  sees  that  it  has  eigenvalues  on  the  imaginary 
axis  along  a  curve  approximately  given  by  a  *  -3e/4. 

We  next  analyze  the  periodic  orbits  of  (3.14)  for  (a,e) 
small.  Any  such  orbit  must  be  close  to  one  of  the  periodic 
orbits  for  (a,e)  =  (0,0).  For  (a,e)  /(0,0)  let  y(t),  p(t) 
be  a  periodic  solution  of  (3.14)  of  period  T,  normalized  so 
that  y(0)  *  0.  Following  the  procedure  in  [3],  one  can  deter¬ 
mine  necessary  and  sufficient  conditions  on  (a,e)  in  order  that 

(3.14)  has  a  periodic  orbit  of  period  T  near  the  orbit  Vp 
defined  by  y(t),  p(t).  This  is  an  application  of  the  Fredholm 
alternative  and  the  method  of  Liapunov-Schmidt  to  obtain  a  bi¬ 
furcation  function  GT(o,e)  with  the  property  that  there  is  such 
a  period  orbit  if  and  only  if  GT(a,e)  =  0.  The  function  GT(a,e) 
has  the  form 

r  T/ 2  2  fT/2  .  2 

(3.15)  GT(a,e)  =  al  p*(t)dt  ♦  ej  p*(t)dt  +  0((|a|  +  |e|r) 

-T/2  -T/2 

as  (a,e)  (0,0). 

The  Implicit  Function  Theorem  implies  there  is  a  unique 
solution  a  *  a*(e,T)  of  (3.15)  for  c  sufficiently  small, 
ot*(0,T)  •  0  and 
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Using  the  first  integral  again,  it  is  not  difficult  to  see  that 
the  period  T  is  a  strictly  decreasing  function  of  p(0)  approach¬ 
ing  2 tr  as  p  -►  1/2.  This  implies  there  is  a  unique  periodic 
orbit  of  period  T  on  each  of  the  curves  defined  by  a  =  a*(e,T) 
above. 

Using  the  differential  equation  and  the  first  integral  V, 
one  can  show  that 


~  f  o2'1!'  J«T  \ 


p (T/2 )  r2 

p(0)  [aTT+r4/3-rZ/4]^2 


dr 


where  aT  *  V(p(0),y(0))  =  p(0)/4  -  p3(0)/12. 

As  T  +  «>,  it  is  not  difficult,  to  see  that  a^T  0, 

6T  -*■  5  as  T  -*■  <o  where  5  ■  -1/2  is  the  number  computed  before 
representing  the  slope  of  the  curve  defining  the  homoclinic  orbit. 
As  T  -*■  2rr ,  aT  -*•  11/96,  6T  -*■  -3/4,  the  slope  of  the  curve 
defining  the  Hopf  bifurcation.  Consequently,  there  is  a  unique 
periodic  orbit  of  (3.14)  between  these  two  extreme  curves. 

We  have,  therefore,  obtained  a  complete  analysis  of  the 
behavior  of  the  solutions  of  (3.14)  for  all  a  6  1R,  e  small. 
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If  we  return  to  the  original  coordinate  system  in  (3.12),  we  can 
state  the  following  result. 

Theorem  3.1,  There  is  a  neighborhood  U  of  (p,y)  *  (0,0)  and 
a  neighborhood  V  of  (X,8)  ■  0  such  that  the  neighborhood  V 
is  divided  into  regions  as  shown  in  Fig.  4  such  that  the  flow  for 
Eg.  (3.14)  in  each  region  is  the  one  depicted  in  Fig.  5.  The 
curves  ri»r2,r3  are  given  approximately  by 


Tj  :  X-v  -$  -B2/2 
r2:  X'v  -B  -  3B2/4 

r3:  X~  -  2$ 


Figure  4. 
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All  curves  in  Fig.  4  correspond  to  saddle-node  bifurcations 
except  T1  and  r2.  Each  point  on  r2  corresponds  to  a  generic 
Hopf  bifurcation  and  each  point  on  corresponds  to  bifurcation 

through  a  homoclinic  orbit. 

We  remark  that  the  same  conclusions  as  in  Theorem  3.1  will 

be  valid  for  (3.11)  for  most  values  of  the  modal  parameters  d,e,f,g. 

2 

Only  the  terms  in  $  of  the  curves  , r2  will  be  changed. 

There  will  be  a  submanifold  of  codimension  2  where  the  complete 
description  of  the  flow  will  require  the  terms  of  order  higher 
than  in  (3.11). 

Theorem  3.1  gives  a  complete  description  of  the  flow  for  the 
autonomous  part  of  Eq.  (3.11)  obtained  by  neglecting  terms  of 
ordeF  >3.  These  higher  order  terms  are  periodic  in  the  independent 
variable.  Thus,  we  have  a  problem  similar  to  the  one  described  in 
the  previous  section.  The  results  there  and  especially  Theorem  2.1 
gives  a  description  of  the  flow  for  the  complete  equations  except 
in  a  neighborhood  of  the  curve  of-  bifurcation  through  a  homo¬ 

clinic  orbit.  The  periodic  perturbation  near  points  on  this  curve 
can  change  the  structure  of  the  flow  in  a  significant  way  as  we 
have  seen  in  Section  2.  Points  on  this  curve  are  analyzed  in  the 
same  manner  as  in  Section  2  in  the  following  way. 

Suppose  we  consider  the  scaled  equations  (3,13).  Applying 
the  method  of  Section  2  to  obtain  the  curves  in  the  a,e  plane 
for  which  there  is  a  homoclinic  orbit,  we  choose  a  point  (oq,Eq). 


corresponding  to  a  point  on  T^,  introduce  a  phase  shift  $ 
along  the  homoclinic  orbit  corresponding  to  (ag,e0).  If 
a  *  dp  +  v,  e  *  £q  +  p,  then  one  can  obtain  a  bifurcation 
function  G(<f>,v,p)  for  C  1R»  (p,v)  close  to  zero.  The 
function  G  has  the  form 

-00  .  QO 

G(<J>,v,u)  =  v  Pq  +  M I  pj  +  e0Mha  E  ( <p )  +  0((|v|  +  |p|)2) 

as  v,p  -*•  0.  If  we  make  the  generic  hypothesis  that  h  (<f>) 

°0  *  e0 

have  absolute  local  maximum  and  minimum  and  no  ther  extreme 

values,  then  there  are  sectors  in  the  (\>,p) -plane  in  which 

there  are  either  no  homoclinic  orbits  or  there  are  homoclinic 

orbits,  the  same  situation  that  occurs  in  Theorem  2.2. 

The  function  h  (<t») ,  in  principle,  can  be  computed. 

0’  0 

One  'should  be  able  to  use  the 

same  argument  as  in  [3]  to  show  that  subharmonic  bifurcations 
also  occur  near  points  on  r^.  Due  to  the  complexity  of  the 
computations,  we  do  not  dwell  on  this  question. 

Equations  (8.8)  are  the  generic  situation  for  two  purely 
imaginary  and  one  zero  eigenvalue.  If  further  symmetries  occur 
in  the  problem,  there  may  be  no  second  order  terms  in  the  normal 
form  for  the  vector  field.  In  this  case,  the  simplest  case  for 
the  approximate  equations  are 
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(3.16) 


P.a  p(X-ap2-by2) 
y  =  y(B+cp2+dy2) 


with  a,b,c,d  fixed  nonzero  constants  and  X,$  small  bifurca¬ 
tion  parameters.  In  the  case  of  a  fourth  order  equation  with 
two  purely  imaginary  roots,  the  same  equations  occur  in  a  natural 
way  coupled  with  two  angle  variables.  The  complete  bifurcation 
diagram  for  this  equation  is  more  difficult  to  obtain  and  will 
require  the  consideration  of  the  fifth  degree  terms. 

The  only  case  that  will  be  discussed  in  detail  is  a  =  d  =  1, 
be  >  1.  We  will  make  some  remarks  on  the  other  cases  later. 

Thus,,  we  consider  the  equations 


(3.17) 


p  *  p(X-p2-by2) 
y  =  y (fJ+cp2+y2) . 


The  results  are  summarized  in  the  following  theorem. 


Theorem  3.2.  There  is  a  neighborhood  U  o£  (p,y)  =  (0,0)  and 
a  neighborhood  V  of  (X,B)  »  (0,0)  •  such  that  the  neighborhood 
V  is  divided  into  regions  as  shown  in  Fig.  6  such  that  the  flow 
for  Eq.  (3.17)  £or_  p  >  0,  y  >  0,  in  each  region  is  the  one 
depicted  in  Fig.  7.  The  curves  rlfr2,r3  are  given  by 
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rx  :  x  ♦  be  *  o,  6  <  o 

r2  :  X (1+c)  +  6 (1+b)  =0,  e  <  0 

Tj  :  Xc  +  B  a  0,  8  £  0 

All  of  the  bifurcations  are  saddle-node  type  except  the  one  that 

takes  place  between  region  4  and  region  5.  On  the  curve 

the  system  is  Hamiltonian  and  the  flow  is  depicted  in  Figure  8. 


Figure  6. 


* 
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Remark.  The  nature  of  the  bifurcation  from  Region  4  to  Region  5 
is  very  complicated  and  not  generic.  To  obtain  a  generic  bifurca¬ 
tion,  one  must  have  some  terms  of  order  higher  than  three.  This 
situation  will  be  discussed  later. 

Proof :  The  lines  p  =  0  and  y  =  0  are  invariant  for  all  X,B- 
In  Regions  1,2, 3, 6, 7,  it  is  easy  to  verify  that  all  equilibrium 
points  lie  on  these  lines.  An  analysis  of  the  linear  variational 
equation  around  these  points  gives  the  flows  depicted  in  these 
regions.  Furthermore,  on  the  curves  r^,r2  there  is  a  saddle-node 
bifurcation  which  gives  rise  to  an  equilibrium  point  (ro»yo^  with 
p0  >  °»  ^0  >  °*  Thus,  if  there  is  to  be  a  periodic  orbit  of  (3.17), 
it  must  occur  for  values  of  A,B  in  the  union  of  the  regions  four 
and'five  and  -T^.  Also,  for  a  periodic  orbit  to  exist,  there  must 
be  some  value  of  a,B  for  which  the  eigenvalues  of  the  linear 
variational  equation  about  (rg,y0)  are  purely  imaginary.  If  A 
is  the  matrix  of  this  equation,  then 

det  A  =  4pgyg (bc-1) 

(3.18)  tr  A  -  -2(pJ-yJ) 

Pg  -  -  (A+Bb)/  (bc-1) ,  y§  =  (B+Ac)/  (bc-1) 

2  2 

Since  be  >  1 ,  pQ  >  0 ,  yQ  >  0 ,  the  eigenvalues  are  purely  imaginary 
if  and  only  if  tr  A  -  0,  that  is, 
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(3.19) 


X(l+c)  +  B(l+b)  *  0 


Equation  (3.19)  is  the  formula  for  The  curves  ,  r2  are 

2  2 

obtained  respectively  when  pjj  ■  0  and  *  0. 

The  local  flow  near  rj,r2  easily  shown  to  be  the  one 
depicted  in  Regions  4  and  5.  This  picture  remains  global  in  both 
regions  4  and  5.  In  fact,  any  periodic  orbit  must  enclose  (pq^q) 
since  there  are  no  other  critical  points.  Such  an  orbit  could 
only  appear  by  the  introduction  of  a  saddle-node  type  bifurcation 
in  region  p  >  0,  y  >  0.  This  is  impossible. 

To  show  there  is  a  first  integral  on  r2  *  {(X,p):X  =  £  0 

mQ  =  (l+b)/(l+c)  >  0,  let  X  *  -mQ& 

(3.20-)  .  P  |  0 1 1/2p,  y  |  B 1 1/2y ,  t  w-  |  8 1  ”*t . 


Since  B  <_  0,  we  obtain  the  equivalent  equations 


(3.21) 


*  2  >  3 

p  =  mQp  -  -bpy  p 

y  *  -y  +  cp2y  +  y3 


If 


.  b+1 

BFT 


v 


c+1 

HcTT 


(3.22) 


V 


then  the  function 


(3.23)  H(u,v)  =  uVll-J  u-v) 

is  a  first  integral  of  Eq.  (3.21).  This  completes  the  proof  of 
the  Theorem. 

Remark.  If  b  >  0 ,  c  >  0 ,  be  <  1 ,  then  the  equilibrium  point 
(p0,y0)  with  p0  >  0,  yQ  >  0  is  always  a  node,  the  curve  Tj 
does  not  occur  and  the  flow  in  the  region  between  and 

is  pictorially  the  same  as  shown  in  Regions  4  and  S  in  Fig.  7 
except  the  interior  point  is  a  stable  hyperbolic  node.  If  b  >  0, 
c  <  0,  then  additional  complications  arise  because  there  is  the 
possibility  of  two  more  equilibrium  points  in  the  region  p  >  0, 
y  >  '6  and  there  can  be  a  saddle-node  bifurcation  in  this  region. 

We  do  not  discuss  this  case  in  detail. 

Let  us  now  discuss  perturbations  6f  Eq.  (3.17),  beginning 
with  autonomous  perturbations.  All  of  the  phase  portraits  in 
Fig.  7  will  be  preserved  under  small  perturbations.  However,  the 
one  in  Fig.  8  will  not.  Therefore  to  obtain  a  more  complete  picture 
of  the  behavior  of  the  equation  under.. perturbation ,  we  consider 
the  effect  of  the  addition  of  higher  order  terms  in  the  same  way 
as  cubic  terms  were  considered  in  Theorem  3.1. 
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If  we  consider  Eq.  (3.16)  as  arising  from  the  problem  in 
1R4  where  the  linear  variational  equation  has  two  pair  of  purely 
imaginary  roots,  then  the  perturbation  terms  must  be  5th  degree. 
Thus,  we  consider 


(3.24) 


p  =  p(X-p2-by2+p4) 
y  =  y(6+cp2+y2). 


We  are  interested  in  the  behavior  of  solutions  near  the  line 
X  =  -m08,  mQ  =  (b+l)/(l+c).  Introduce  the  scaling 


6|1/2p,  y 


|e|1/2y,  t  |  8 1  *At ,  X  H-  -m0B  ♦  aB 


-1. 


to  obtain 
♦ 

(3.25) 


P  *  p(m0+a-p2-by2+Bp4) 
y  *  y(-l+cp2+y2).  . 


One  can  proceed  as  in  the  proof  of  Theorem  3.1  to  obtain  the 
curve  in  the  (X, 8) -plane  along  which  Eq.  (3.24)  has  a  homoclinic 
orbit.  In  terms  of  the  scaled  variables,  we  observe  that  such  a 
solution  can  exist  only  in  a  neighborhood  of  a  *  0,  6  *  0.  We 
then  obtain  a. bifurcation  function  G(o,8)  for  a  homoclinic  orbit 
and  observe  that  it  has  the  form 
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G(a,0)  -  of  p2(t)dt  ♦  ef  Pg  (t)dt  +  0((|a|  +  |6|)2) 

J  .  OO  v  -00 

as  a,0  •+■  0,  where  (Pq (t)  ,yQ (t) }  describe  the  homoclinic  orbit 
for  a  *  0,  0*0.  There  is  no  term  involving  pg(t)  because 
for  a*0,  6*0,  the  system  is  Hamiltonian.  The  equation 
G(o,6)  *  0  has  a  unique  solution  a  *  a*(B)  in  a  neighborhood 
of  (a, 0)  *  {3,0),  a*(0)  *  0.  This  gives  a  curve 

t'2  -  { (X, B)  :  A  -  -  m0B  ♦  o*(0)0) 

f 

in  the  (A,0)-plane  such  that  Eq.  (3.4)  for  (A,0)  €  r2  has  a 
homoclinic  orbit. 

M 

_0ne  can  also  obtain  a  curve  ?2  in  the  ( X , 6 ) -plane 
T2  *  { ( X , B)  :  X  *  -  m08  +  o**(6)6) 

t» 

where  o**(0)  *  0,  such  that  for  any  ( X , B)  €  r2,  the  linear 
variational  equation  for  the  solution  (PqjYo)*  Po  ?  >  ®» 

of  (3.4)  has  both  eigenvalues  purely  imaginary.  Proceeding  as 
before,  one  obtains  the  existence  of  a  unique  periodic  orbit  in 
the  Region  4*  between  r2  and  r2.  The  flow  in  Region  4'  is 
shown  in  Fig.  9.  The  complete  bifurcation  diagram  is  shown  in 
Fig.  10. 


It  remains  to  discuss  perturbations  of  (3. 24) when  the 
higher  order  terms  contain  angle  variables.  If  there  is  only 
one  angle  variable  (the  case  of  two  purely  imaginary  and  one  zero 
root)  everything  follows  exactly  as  for  Eq.  (3.14).  If  there  are 
two  angle  variables,  the  situation  is  more  complicated.  The 
general  theory  of  integral  manifolds  implies  that  every  hyperbolic 
equilibrium  point  of  the  unperturbed  equation  becomes  a  hyperbolic 
torus  of  dimension  2  and  every  hyperbolic  periodic  orbit  of  the 
unperturbed  equation  becomes  a  hyperbolic  three  dimensional  torus 
as  one  crosses  ^  (the  Hopf  bifurcation).  The  analysis  in  Section 
2  shows  that  the  two  dimensional  torus  bifurcates  to  give  a  three 
dimensional  torus.  None  of  the  other  cases  are  well  understood. 
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